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Abstract. In this paper we consider the classical initial value problem for 
the bosonic membrane in light cone gauge. A Hamiltonian reduction gives a 
system with one constraint, the area preserving constraint. The Hamiltonian 
evolution equations corresponding to this system, however, fail to be hyper- 
bolic. Making use of the area preserving constraint, an equivalent system of 
evolution equations is found, which is hyperbolic and has a well-posed initial 
value problem. We are thus able to solve the initial value problem for the 
Hamiltonian evolution equations by means of this equivalent system. We fur- 
thermore obtain a blowup criterion for the membrane evolution equations, and 
show, making use of the constraint, that one may achieve improved regularity 
estimates. 



1. Introduction 

The initial value problem for the classical evolution of a physical system seeks 
to characterize critical points of the action functional associated to the problem 
in terms of an appropriate set of initial data. Once the well-posedness of the 
classical field equations has been shown, the initial data not only determines the 
classical motion for some time interval (0,T), where T is the time of existence 
corresponding to the data, but also the Hilbert space of wave functions of the 
corresponding quantum-mechanical system. For physical systems without gauge 
symmetries described by conjugate pairs [x,p) satisfying Hamilton's equations, the 
initial data is given directly by specifying the conjugate pair (ccqjPo) at an initial 
time. The wave functions, in the Schrodinger picture, are precisely the space of 
functions (/?(xo) with ip E LF' . 

In the presence of gauge symmetries the initial data is restricted by constraints 
and gauge- fixing conditions. One may solve these restrictions at the classical level 
in terms of conjugate pairs and then determine the wave function as before, or 
one may consider general wave functions lp{xq), Lp £ L'^, and restrict them at the 
quantum level to a subspace !K C ; the domain of the quantum operators should 
then be dense in "K. In this latter case, one may then extend the phase space in 
order to realize the BRST symmetry of the quantum system. A key point in both 
procedures is to determine both the classical and quantum restrictions associated 
to the gauge symmetries. 

The formulation of the initial value problem for gauge theories, including Ein- 
stein gravity, electromagnetism, Yang-Mills, string, and membrane theories, as well 
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as their supersymmetric extensions, is formally solved by the approach of Dirac 
[l^Fl It determines in a constructive and systematic manner, the constraints asso- 
ciated to the gauge theory. Moreover the procedure ensures that the constraints are 
preserved in time by the Hamiltonian flow, with the Lagrange multipliers associated 
second-class constraints determined by the conservation procedure and Lagrange 
multipliers associated to first-class constraints remaining as gauge-dependent vari- 
ables. The resulting Hamiltonian formulation of the gauge theory ensures that if 
the constraints are satisfied initially, then they are also satisfied on (0, T), for some 
time of existence T. The first-class constraints close as an algebra under a Poisson 
bracket constructed from the symplectic structure of the Hamiltonian formulation. 
The generators of this algebra realize the gauge symmetry of the action functional, 
and consequently of the field equations themselves. (For a detailed exposition of 
this systematic approach in the case the Einstein gravity, including a discussion of 
first- and second-class constraints, see (30|.) 

A well-known gauge choice of the above mentioned field equations is the light 
cone gauge (LCG). In this gauge, one can solve the constraints, together with the 
gauge fixing conditions, in terms of unconstrained physical degrees of freedom. 
This property of the LCG becomes very useful when proving relevant properties of 
the corresponding quantum theories, one example being arguments concerning the 
unitarity of iS-matrix in superstring theory. A treatment of Einstein gravitation has 



been considered in [43|, [3l|, In [J] the positivity of the energy and couphng 
to gauge fields was analyzed. For a review of non-covariant gauges in gauge field 
theories, including the light cone gauge, see 0|. 

The LCG has proven particularl y fr uitful in classical and quantum analyses of 



the D — 11 supermembrane theory , [27[ , [16| , which is also a a relevant ingredi 



ent of M-theory. The supermembrane theory has first-class constraints associated 
to the generators of diffeomorphisms of the world volume and the local fermionic 
symmetry, known as K-symmetry, as well as second-class fermionic constraints. The 
complete set of constraints is very difficult to treat in a covariant formulation. How- 
ever, the LCG allows an explicit solution of all first- and second-class constraints. 
It is also convenient in the case of both membrane and supermembrane theory 
to fix all symmetries up to area-preserving diffeomorphisms, which in this case 
are in fact symplectomorphisms (with respect to a symplectic structure defined as 
part of the gauge choice) 3 The resulting Hamiltonian, in LCG with these resid- 
ual symmetries, may be analyzed without difficulty. The constraint associated to 
the area-preserving diffeomorphisms may be interpreted as a symplectic generaliza- 
tion of the Gauss law, with nonlinear terms a la Yang-Mills, but arising from the 
symplectic bracket rather than the bracket of the SU{n) Lie algebra ^38;]. 

Once the LCG is implemented in a Hamiltonian formulation of the supermem- 
brane, and the first- and second-class constraints are solved, one obtains directly a 
canonical Hamiltonian reduction of the original formulation. That is, the elimina- 
tion of phase-space variables occurs in canonical-conjugate pairs. The Hamiltonian 



For the formal computations to be properly defined, the well-posedness of the Hamiltonian equa- 
tions must be established. 

■^In even dimensions different from two, these symplectomorphisms are volume-preserving, but not 
all volume-preserving diffeomorphisms are symplectomorphisms. 
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in the LCG is polynomial in the remaining variables and their derivatives and, be- 
cause it is formulated in terms of physical degrees of freedom (the residual gauge 
symmetry may be fixed in a convenient manner), it is essentially the same Hamil- 
tonian which appears in the path-integral formulation of the quantum theory. This 
means that properties of the same potential will determine both classical and quan- 
tum aspects of the theory. The classical and quantum stability properties of the 
membrane or supermembrane theories are determined by the nonlinear dependence 
of the potential, along the configurations for which the potential becomes zero. 

These properties have been analyzed by considering a SU (n) regularization of 



the membrane or supermembrane theory 27[ , 17| , [l8l | . This regularization is itself 



an interesting physical model, and was the starting point in the introduction of the 
matrix model. The S'C/(n)-regularized supermembrane is a maximally supersym- 
metric Yang-Mills theory in 1 + dimensions [3] , the classical field equations being 
ordinary differential equations and the corresponding quantum problem finite di- 
mensional. The regularized membrane potential has "valleys" extending to infinity; 
the value of the potential at the bottom of these valleys is zero. Thus the static 
solutions of the equations of motion with zero value of the potential are unstable. 
Nonetheless, the quantum-mechanical Hamiltonian of the membrane theory has a 
discrete spectrum, due to the structure of the valleys. In dimensions greater than 
one, the discreteness of the spectrum of the regularized Hamiltonian, a Schrodinger 
operator, is determined by the behavior of the mean value of the potential in the 
sense of Molcanov [4l| , [s^l . This mean value tends to infinity as one moves out- 
ward along the valleys in configuration space, thus ensuring that the operator has 



discrete spectrum |47[,[36|, [25|. In the supermembrane case, the potential becomes 
unbounded from below, due to Fermionic contributions, in a manner which renders 
the spectrum continuous ^17J . However, the supermembrane with central charges 
generated by the wrapping of the supermembrane on a compact sector of the tar- 
get space has discrete spectrum ^lij . In fact the topological condition ensuring the 
nontrivial wrapping, which does not modify the number of local degrees of freedom, 
eliminates the nontrivial configurations with zero potential. 

Besides this interesting relation between classical and quantum stability prop- 
erties of the membrane and supermembrane theories, there are other aspects of 
classical supermembrane theory which reproduce quantum a' effects of string the- 
ory. (That is, perturbative quantum effects in string theory where the perturbative 
parameter is the inverse of the string tension.) The closure of the k symmetry in 
supermembrane theories with a general background metric on the target space is 
only possible provided the background satisfies the D = 11 supergravity equations 
|12| . The analogous result for superstring theory arises only when a' quantum 
effects are taken into account. Furthermore, the HA _D-brane action in D = 10 
may be obtained from a duality transformation from the D = 11 supermembrane 
compactified on a circle. The same _D-brane action arises from Dirichlet strings 



only when quantum effects are considered [44| . [4 



In light of these considerations, it is natural to analyze the classical initial- value 
problem of the membrane and supermembrane in the LCG. This not only provides 
a foundation for the study of the quantum mechanical systems corresponding to 
the membrane and supermembrane in LCG, by putting the classical theory on a 
firm basis but, in establishing a criterion for continuing the classical solution in 
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time, is an important first step towards the identification and characterization of 
whatever singularities may develop. This is interesting not only from a classical 
perspective, but also from a quantum mechanical one, as singularities are expected 
to play an important role in the quantum theory It may also provide a 

framework to analyze the large n limit of the regularized theories and the related 
stability problems, one may hope to extrapolate consequences for the quantum 
stability problem. 

From a mathematical perspective, the problem is interesting as the membrane 
equations, which correspond to the supermembrane field equations after the spinor 
dependence (i.e., the fermionic sector) has been anhilated, are one case of an im- 
portant class of geometric wave equations. Geometrically, membranes are timelike 
submanifolds with vanishing mean curvature; the equation governing this condition 
is the Lorentzian analogue of the minimal submanfold equations (much as wave 
maps are the Lorentzian analogue of harmonic maps) . Attempts to approach the 
problem of existence of such submanifolds by applying techniques from the theory 
of differential equations are complicated by the inherent diffeomorphism invariance 
of the problem; in an arbitrary coordinate system the equations are not strictly 
hyperbolic (i.e., wave equations). As in the case of the mathematical study of the 
Einstein field equations, this difficulty can be overcome by choosing a gauge, which 
eliminates (or at least reduces) the diffeomorphism freedom and yields a system of 
equations to which PDE theory can be applied. 

In fact, there are actually two levels at which the mathematical problem of local 
existence can be posed. As the equations governing the embedding are wave-type 
equations, one expects to pose an initial value problem. The first is at the level 
of geometry: Given an initial spacelike submanifold and timelike vectorfield on the 
submanifold, can one extend the submanifold in the direction of the vectorfield such 
that the mean curvature of the extension vanishes? The second is at the level of of 
embedding functions: Given an embedding function for an initial spacelike slice, an 
initial 'velocity' for that function, and a gauge condition, can one find an embedding 
function satisfying the relevant PDE (as expressed under the gauge condition)? In 
either case, one would also like to show Cauchy stability as well: that not only do 
local solutions exist, but that they are unique and depend continuously upon the 
given data. 

The problem of local well-posedness at the level of geometry has been recently 
addressed in a very general setting in [i^ . At the level of the PDE, that work makes 
use of a variation of the harmonic gauge condition, which was first applied to the 
membrane problem in a Minkowski ambient spacetime by and has been used 
extensively in the study of the initial value problem for the Einstein equations ([13], 
see also |22|). These works provide a very satisfactory resolution to the geometric 
local-existence problem, as well as a rather complete solution to the local existence 
problem for the PDE in harmonic coordinates. It leaves open, however, the problem 
of local existence of embedding functions satisfying other gauge conditions which 
may be better-suited for addressing questions of the lifespan of solutions (and/or 
singularity formation of solutions) or questions arising when considering aspects of 
the quantum problem. 
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In fact, relatively little is known concerning the existence of solutions to the PDE 
when reduced by even the simplest gauge conditions. A Hamiltonian reduction 
under the partial gauge condition that the time coordinate of the submanifold 
coincide with the Minkowski time coordinate was considered in 42| (see also 26]). 
Under this choice of foliation, the equations for codimension-1 membranes reduce 
to a first-order system. A large number of examples of solutions satisfying this 
gauge condition have been constructed; see for example [28| . 

In this work we address the well-posedness of the membrane field equations in the 
LCG. In particular, we show that for any non-degenerate initial data satisfying the 
constraints there exists a time interval (0, T), with T depending on the initial data, 
and a unique solution to equations of motion of the reduced (in LCG) Hamiltonian 
corresponding to the initial data. Furthermore, both T and the solution depend 
continuously on the initial data in a suitable topology. The result is obtained by 
application of the theory of hyperbolic partial differential equations. 

The system of field equations obtained by taking variations of the membrane 
action in light cone gauge is a fully non-linear system which is not, however, hy- 
perbolic. We therefore introduce a modified system which is quasi-linear and hy- 
perbolic, and has the property that solutions of this modified system satisfy the 
constraints and field equations if these are satisfied initially. It would be interesting 
if this modified system has a (super-) membrane action associated to it, as it has 
additional constraints which are preserved by evolution and which should appear 
in any quantum formulation of such a theory. The existence of such preserved con- 
straints is an indication that there may be a gauge theory, with gauge symmetries 
generated by both constraints, which under some appropriate gauge fixing reduces 
to the modified system we employ. 

The modified system is obtained by differentiating the LCG field equations with 
respect to time, and eliminating terms which vanish due to the constraints. Similar 
techniques have been used to extract hyperbolic systems for the Einstein and Yang- 
Mills equations without gauge fixing, see [l| and references therein. We are able to 
show the existence of solutions to the modified equation using a standard argument 
based on energy estimates. The structure of the modified system is such that 
the normal procedure for obtaining energy estimates, commuting spatial derivative 
operators through the equation, leads to a loss (see commutator estimate (jS.llcp ). 
We recover this loss by means of an elliptic estimate and by commuting a time 
derivative through the equation. The result of our method is a local existence result 
which requires slightly more regularity than expected for arguments based on the 
Sobolev embedding. We are nevertheless able, by making use of the constraint, to 
obtain an improved energy estimate for solutions the original equation. With this 
improved energy estimate we are able to show that the time of existence depends 
on the "classically-expected" norm of the initial data. 

In this paper we work in terms of integer order Sobolev spaces only. By making 
use of more sophisticated techniques, the results presented here can be improved 
as far as the regularity requirements are concerned. The algebraic structure of the 
reduced field equation, makes it interesting to ask for the optimal well-posedness 
result from the point of view of the regularity of initial data. The matrix analog of 
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the system gives an ODE analogof the membrane system, which has been exten- 
sively studied (see for example [3] , [1^ , ) • It is interesting to consider this from 
an analytical point of view as a consistent truncation of the system, and to make 
use of related ideas to study the local well-posedness of the system for rough initial 
data. 

The organization of the remainder of this paper is as follows. In the next section 
we outline our notational conventions, and list a number of functions spaces and 
related estimates appearing in the local existence proof. In subsequent section [31 
we introduce the Lagrangian formulation of the membrane problem and perform 
a Dirac-style canonical analysis of the membrane problem, deriving the reduced 
equations of motion under the light cone gauge condition. In section |4] we give a 
treatment of the initial value problem by means of the modified system described 
above. The modified system is hyperbolic and well-posedness follows along essen- 
tiall y st andard lines, with the additional difficulty that the system is fully nonlinear. 



See [32|, [35| for treatments of related problems. For completeness, we give a self- 
contained proof in section [5] Finally in section [6] we derive the improved energy 
estimate, which gives us the improved estimate for the time of existence. 



2. Preliminaries 

We consider 3-dimensional submanifolds M of _D-dimensional Minkowski space 
with M = E X E and S some compact 2-manifold. 

We make use of a number of index sets: Greek indices /i, z^, . . . , ranging over 
0, — 1, refer to Cartesian coordinates in Minkowski space. Middle Latin 

indices m, n, . . . , ranging 1, . . . , D — 2, refer to Cartesian coordinates on Euclidean 
space R-^^^. Lower-case early Latin indices a, 6, ... , take values 1, 2 and refer to 
coordinates on S. Upper-case Latin indices A, B, . . . , take values 0, 1, 2 and refer 
to coordinates on M . We now describe the various coordinates in more detail and 
indicate which metrics are used to raise/lower each set of indices. In all cases we 
sum over repeated indices. 

In Cartesian coordinates {x^) the Minkowski metric r;^^ is given by 

T]f,^dx''dx'' = -{dx")"^ + [dx^f + ■■■ + {dx^-^f. (2.1) 

Greek are raised/lowered using ry^^ so that x^ = rj^j^^x'^ . 

We also make use of null coordinates x^ = {x'^ ± x^^^) for Minkowski space, 

denoting by (a;™), the remaining coordinates (a;^, . . . , in Euclidean space 

y^f^iY^ respect to these coordinates the Minkowski metric is given by = 
rj — 0, 7/-I = rj y- — —1, and rjmn = 5mn is the Euclidean metric. 

Let T be some global coordinate whose level sets S,- foliate M . When there is no 
confusion, we denote Y^r by S. On each E we use (r-independcnt) local coordinates 
((T°) with a = 1,2. Together (r, cr") = (^^) give coordinates on M. 

Any (embedded) submanifold M C MP is determined by the function x = {x^) : 
M R-^ which induces a metric gAB on M which in local coordinates is given by 
gAB — Vtiyf^AX^dBx"] here — |fx. In what follows we restrict attention to 
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the case where the induced metric gAB is Lorentzian with dr a timehke direction. 
The metric qab induces a volume form fj,g, given in coordinates by /ig = rf^*^ A 
d^^ A with y/\g\ = ^ — det \gAB[ , used below to define the Lagrangian action 

for the membrane system. 

We denote by jab the metric on S induced by gAB', we require jab to be Rie- 
mannian. Note that 

lab = VniydaX'^dbX" = -{daX~^dbX~ + daX~ dbX~^) + daX™dbXm- (2.2) 

Where we desire to indicate explicitly the x-dependence of jab we write j{x)ab- The 
metric jab gives rise to a volume form /i-y = ^fj d<j^ Ada'^; here 7 is the determinant 
of Jab- By the usual formula for the inverse of a matrix, the inverse 7"^ of jab can 
be expressed 

^ab ^ 1 gacgbd (2.3) 

7 

Here is the anti-symmetric symbol with two indices. (Explicitly, e^^= — ^^^ = 
1, g^^=g^2= 0.) Using the anti-symmetric symbol, the determinant 7 can be 
written as 

7=1 (2.4) 

2.1. Symplectic structure. The light cone gauge condition requires the choice of 

a fixed (r-independent) area form ^Jwda^ A da"^ on S. We presume ^/wda^ A da'^ 
to be the area element arising from some fixed background metric Wab on E; thus 
^/w = Vdet Wab- We may without loss of generality assume that Wab is real analytic. 

The area form gives rise to a symplectic structure on S. In local coordinates, 
the Poisson bracket associated to this symplectic structure is given by 

{f,g} = ^dafdbg (2.5) 
w 



where f,g are functions on E. Recall that the Poisson bracket is bilinear, skew 
{{fi9} — ^■ud satisfies the Jacobi identity 

= {/, {5, h}} + {g, {h, /}} + {/i, {/, g}}. (2.6) 

In what follows we make use of the area form ^/w da^ A da"^ when integrating on E, 
but denote the form simply by ^/w (i.e. the expression da^ Ada^ is to be implicitly 
understood). Note that Stokes' theorem implies 

/ {f,g}h^f^=- [ f{h,g}V^. (2.7) 

2.2. Derivatives and function spaces. Here we outline our notational conven- 
tions regarding derivatives and introduce some function spaces used in the local 
existence results below. These spaces are defined with respect to a fixed atlas of 
coordinate charts on E. 

For function v defined on [0, T] x E, we denote by Dv any coordinate derivative 
daV. By D'"v, with I a non-negative integer, we mean an arbitrary combination 
of I coordinate derivatives of v. We furthermore denote by dv any of daV or drV. 
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In a slight abuse of notation, in the presence of norms we implicitly sum over all 
coordinate derivatives, so that 

\dvf = \drvf + \Dvf; (2.8) 

\D''v\, \Ddv\, etc. are defined analogously. 

The following is an overview of the norms and function spaces used. 

Sobolev spaces: Denote by il' the Sobolev space of functions on S whose 
derivatives (in a fixed coordinate atlas) of up to order I are square-integrable. 
Thus 

Mh= f kpy^, \\v\\l. = \\v\\h+ [ \Dv\'V^, etc. (2.9) 

L°° spaces: The following norms control the (essential) supremum of func- 
tions: 

||i;||loo = esssup|u|. (2.10) 

E 

When V S C°(S) one can replace esssupj^ with maxs- Furthermore define 

\\v\\w^.^ = \\v\\L^ + \\Dvh^. (2.11) 

Curves in function spaces: The set of maps w : [0,T] which are r 

times differentiable with respect to r is denoted C"'([0, T]; iJ') and given 
the norm 

r 

\\w\\cr-m,T];H^) = sup V||a>||H!. 

For functions w G ^([O, T]; ff'), we use a subscript to denote restriction 
to T = 0. It is important that this restriction is made after any differenti- 
ation, so that (for example) |9woP = |9rU)(0)p + |Dw(0)p. 
Spacetime norms: We make use of two particular norms for curves of 
functions. For 

2 

V e&j,:=P\C\[Q,T\-H^-') (2.12) 

i=l 

we make use of the norms 

NI?=Ell^>ll?.-. (2.13) 

and 

|||^;|||i,. = sup||i;(-)|b. (2.14) 

[0,r] 

Note that is equivalent to + ||99T-?;||i/i-2. We also use the norms 

= Mil + \\drx\\% (2.15) 

and 

{{{x)))i,r = supiixh (2.16) 

[0,r] 

forxGC\[0,T];H^). 
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2.3. Results from analysis. We make use of the following estimates, which can 
be proven using classical methods of calculus; see, for example, Chapter 13 §3 of 
(4^ . Note that while versions of these estimates hold in all dimensions, as presented 
here the estimates are dependent on the dimension of S being 2. Here, and in the 
application of these estimates below, C is a constant independent of the function(s) 
being estimated (and unless otherwise specified depends only on E, our coordinate 
charts, \/w , and the number of derivatives being estimated). 

Sobolev inequality: For / > 1 and v G H' we have w € C° and 

||f||L~ <C||t;||ff,. (2.17) 
Product estimate: For v,w £ n we have 

\\vw\\hi <C{\\v\\l^\\w\\hi+\\w\\l^\\v\\hi) (2.18) 

Note that this estimate, together with (|2.17p implies that \\xy\\i < C|la;||;||y||/ 
provided Z > 1. 

Elliptic regularity: For uniformly elliptic operator L = da[a°'^db{-)] + b°^da 
with a G VF1^°°, 6° g L°° , there exists a constant C, depending on the norm 
of the coefficients and the ellipticity constant, such that 

\\v\\h-2<C{\\v\\l2 + \\Lv\\l^) (2.19) 

Gagliardo-Nirenberg-Moser estimate: For I < k we have 

WD'vWr^..,, < C\\v\\]-J\\D'v\\l, (2.20) 

An easy consequence of Holder's inequality and (|2.20p is the following prod- 
uct estimate for u,v G 

\\{Du){Dv)\\l2 < C\\u\\hA\v\\h^- (2.21) 

Commutator estimate: The estimate (j2.20p also implies the following com- 
mutator estimate 

\\[D\v]w\\hi <C {\\w\\LAAm+^ +\\Dv\\l^\\w\\hi+.-^) ■ (2.22) 

3. Canonical analysis and the light cone gauge 

We now perform an ADM-style canonical analysis of the membrane system. 
Beginning with a Lagrangian action integral, we give a brief treatment for the 
membrane system in general before analyzing in detail the reduction under the 
light cone gauge condition. 

3.1. Lagrangian formulation. We seek an embedding 

X = [x^') : M^+'^ , (3.1) 
which is critical with respect to the action given by the induced volume element 

S = - / (3.2) 
JM 
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In string theory (13. 2|) corresponds to the Nambu-Goto action. An equivalent for- 
mulation, expressed in terms of a polynomial Lagrangian density and corresponding 
to the Polyakov action of string theory, is given by 



^ J M 



(3.3) 



where the inverse metric g"^^ is treated as an independent variable. Both (|3.2p and 
(13.31) lead to the same Hamiltionian formulation. 



Critical points of the functional (13. 2|) give rise to submanifolds M C MP with 
vanishing mean curvature; the Euler-Lagrange equations for (j3.2p are 



which can be seen by computing 



0, ^ = 0, ...D-l, 



\9\ 



1 



S9AB = Vlglg^^dAxf'dBiSa 



(3.4) 



(3.5) 



When \/]g\ ^ 0, the system (|3.4I) can be expressed (in any coordinates) as 

(5^. - g^'^'dcx^dDX,) g^^'dAdBx" = 0. (3.6) 



One can interpret l|3.4p as evolution equations for x{t) : E — > R-^. As the 
induced metric gAB is Lorentzian, one expects the system to be hyperbolic (i.e., a 
wave-type equation) and thus to be able to pose the following initial value problem: 
For initial embedding xq, and initial velocity uq does there exists an interval [0,T] 
on which x{t) is defined and satisfies a;(0) = xq, dTx{Q) — uq, and (|3.4p ? 

As expressed above, the evolution equations for x'^ are not strictly hyperbolic; 
this is due to the diffeomorphism invariance of ()3.2|) under re-parameterizations of 
M. (Note that it is also invariant under choice of coordinates in the Minkowski 
spacetime, but we have fixed these degrees of freedom.) Thus we turn to the issue 
of gauge choice by performing a canonical analysis. 



3.2. Canonical (Hamiltonian) analysis. The Hamiltonian approach plays a 
fundamental role in the classical and quantum analysis of field theories. The for- 
mulation for gauge theories briefly described below was developed by Dirac in [l^ . 
For a detailed description of this approach applied to classical covariant field theo- 
ries see [s^ ; for a descrigtion of path- integral quantum analysis in the presence of 



general constraints see [45 1 



The starting point, as proposed by Dirac [l9|, is an action integral expressed in 
terms of a Lagrangian density defined on some time- foliated manifold M x E. 
The density L depends on some number of independent fields (f> and their spatial 
derivatives up to order k, D(j), . . . , D^4>, as well as on the derivatives of the first 
time derivative of </>, 9t(/), Ddtcj), . . . ^D^dtct)- 

One then introduces the Hamiltonian density via a Legendre transformation 

'K= ( {TTdt(t)-L) (3.7) 
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where the conjugate momenta tt associated to 4>, defined as 

5L 



appear in the functional derivative of the Lagrangian action L — /jjxs'^ with 
respect to independent variations 54) and 5{dt(t>): 

The canonical variables (p, tt take values in an infinite-dimensional manifold 
referred to as the phase space of the theory, which is equipped with a Poisson 
structure given (as a density) by 

\uv]'p = ——-—— (3 10) 

5(j) dir 5(f) 5tt 

The symplectic structure determined by the Poisson bracket plays a fundamental 
role in both the classical analysis of the field theory, as well as in the canonical 
quantization of the theory It is also the main algebraic structure in the deformation 
quantization approach ^l3|, [llj, [s^. 

In gauge theories, dt(t> cannot be expressed in terms of unconstrained momenta 
as the Hessian of iL with respect to dt(j) becomes singular; thus there are constraints 
on the phase space for (</>, tt). Further constraints, or restrictions on the Lagrange 
multipliers associated to these constraints, may arise upon imposing the require- 
ment that the vanishing of the constraints be preserved by evolution under the 
Hamiltonian flow (see below). These constraints are of two types. First-class con- 
straints are those which commute, on the constraint submanifold (defined to be 
the submanifold where all first and second class constraints are satisfied), with all 
other constraints; the Lagrange multipliers associated to these first-class constraints 
are gauge-dependent fields and remain undetermined during the canonical analysis. 
Second-class constraints do not commute (on the constraint submanifold) with all 
other constraints and the associated Lagrange multipliers are determined by the 
condition that the constraints be preserved. 

Notice that there is an implicit assumption in the Dirac approach concerning the 
structure of the constraints: they must be regular. A point xq is a regular point of 
(j) : X Y a 0'(xo) is onto. The constraint </) = is regular when each point of the 
constraint submanifold {x : ^{x) = 0} is a regular point of 0. Irregular constraints 
must be treated in a separate manner as the usual theory for Lagrange multipliers 
assumes regular constraints. 

Once the Hamiltonian H = J^'K and the constraints have been determined, one 
can reformulate the action integral in terms of the canonical fields 

§[0,7r]=/ {7Tdt4 - ^ - KC'') , (3.11) 

where C° define the constraints on the phase space and Xa are the associated 
Lagrange multipliers. The canonical Hamiltonian density defined by 

H,= I 'K,= ( 'K + KC (3.12) 
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determines the evolution of the canonical fields via the Hamiltonian field equations 

dt(t> = [0, ffj-p = 

which are obtained by varying p. lip , and are equivalent to the Euler-Lagrange 
field equations associated to the Lagrangian L. Note that for a general quantity 
F = f{(j), TT, t) one has 

dtF=[f,H,]^ + dtf (3.14) 

along the Hamiltonian flow. 

The action p. lip is also the starting point for the Feynman path-integral formu- 
lation of quantum field theory. Under some assumptions on the dependence of !K 
on TT, the path integral defined from §[(/), tt] is formally equivalent to the one defined 
from the Lagrangian action integral. 

The action integral S[0, tt] is invariant under the gauge transformations generated 
by the first class constraints, this can be easily seen by noting that the procedure 
above ensures that any first class constraint C has [C, iJ]-p = on the constraint 
submanifold. (An interesting feature of diffeomorphism-invariant gauge theories is 
that 'K — Q, i. e. !Kc — AqC".) This gauge invariance leads to degeneracies for the 
classical field equations p.l3p ; thus one typically performs a (partial) gauge fixing 
before proceeding to analyze the equations of motion. There is a general method 
for introducing gauge- fixing terms, and the corresponding Faddeev-Popov terms 
in the path-integral formulation, without solving the constraints. The resulting 
effective action becomes BRST-invariant; it may be obtained from the Hamiltonian 
formulation following 24|, or from the Lagrangian formulation following [§| 



Both approaches present difficulties in the presence of complicated second class 
constraints, such as those arising in supermembrane theory. In that case, it is most 
convenient to explicitly solve the constraints at the classical level (in the light cone 
gauge) and then proceed to quantize the theory. Thus wc follow this latter approach 
in this paper. 



3.3. Canonical analysis for membranes. We start from the Lagrangian action 
p.3p . where and g"^^ arc independent fields. Alternatively, one may start from 
p.2p : the resulting Hamiltonians are exactly the same. 

We perform the usual ADM decomposition of g with respect to the foliation S,-, 

denoting by jab the metric on S, by = ^ ■\/— 5°" the lapse, and by A° — 

N^g°°' the shift vector. We raise and lower a, b with ■^ab, 7"''; thus Na = jabN'^, etc. 
To be explicit, the metric gAB and its inverse g"^^ are given by 

, , / -A2 + N^N, \ ^ AB,( -N~^ N-^N^ \ 

mB\ lab ) N-^N^ 7°^ - N'^N^N' J ■ 

(3.15) 
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Treating x'^ as canonical variables, we see that the conjugate momenta to are 
given by 

= ^ (drX^ - N^daX^,) . (3.16) 

Instead of introducing the conjugate momenta to N , Na, and jabi it is convenient 
to treat them as auxihary fields. Furthermore, jab may be eliminated as an inde- 
pendent field and expressed in terms of the using jab — daX^dbX^^. 

We are able to solve (|3.16l) for drX^^ in terms of 

drx^" = — + N^dax^". (3.17) 
The canonical Hamiltonian density is therefore given by 

^^(p'+7) (3.18) 

Here = p^p'' . 

There are two constraints 

<i> = i(p2 + 7) =0 (3.19) 
*a = Pp5,x^ = (3.20) 

with Lagrange multipliers 

A = — and A'^ = 7V^ (3.21) 
The action can now be written 

p^^rx^' - f (A$ + A'^*,) . (3.22) 

M J M 

Note that <i>a is given independently of the target (Minkowski) metric, while $ 
depends on the ambient metric 7]. 

The constraints $, ^a are first-class constraints in the sense of Dirac. The quan- 
tities $, $a are the generators of time and spatial diffeomorphisms, respectively. 

One can compare the structure of this Hamiltonian to that arising in the theory 
of general relativity (see, for example [13]); the Hamiltonian has the same linear 
structure (i.e., it is linear in $, $a) since both are invariant under diffeomorphisms. 
We also point out that $ has the same quadratic dependence on the momentum as 
the corresponding time generator in general relativity. 



3.4. Light cone gauge. We consider a (partial) gauge fixing of the above action: 
the light cone gauge. It has the property that the gauge fixing procedure gives rise 
to a canonical reduction of the above action, and is also the only known gauge 
where the K-symmetry constraints of the supermembrane can be explicitly solved. 

In order to specify the light cone gauge, we make use of the null coordinates 
(a;^, x~ , X™) in Minkowski space and also the (r-independent) volume form ^Jw on 
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The hght cone gauge condition for the membrane Hamiltonian is determined by 
taking 

a;+ = -p° r and p- = p°_y/w , (3.23) 

where p° is some constantH As is made evident befow, this is only a partial gauge 
fixing, the resulting system being invariant under diffeomorphisms which are area- 
preserving with respect to -^/w . Note that under this gauge choice, the metric 
7afc = daX^dhXm] i-G., it does not depend on derivatives of x^. 

We now proceed to construct the (partial) gauge-fixed Hamiltonian in light cone 
gauge. The constraint $ = may be solved algebraically for p+ in terms of a;™ and 

Pm 

P+ = -^\{PraP'^+l), (3.24) 

C\/w 2 



while the constraint $a = determines x (in terms of and pm) via the relation 

daX- = ^^PrndaX"", (3.25) 

Cy/W 

provided the integrability condition 

^dx"' = (3.26) 

c 

holds for all closed curves C in E. The conjugate pairs x~^,p^ and are thus 
eliminated provided this condition holds. 

The Poisson bracket analysis of (|3.26p shows that it is a first-class constraint 
generating area-preserving diffeomorphisms and is equivalent to the local constraint 

^™ x"l=0 (3.27) 



/w 

in combination with the global constraint 



w 



da;" = (3.28) 



on some basis {6^} of the homology of S. Together, (|3.27p - (|3.28p generate area- 
preserving diffeomorphisms homotopic to the identity. The distinction between 
them is the following. The left side of (|3.27p generates area-preserving diffeomor- 
phisms under the Poisson bracket, with infinitesimal parameter e, a time-dependent 
single- valued function on S, i.e. 

5,a;™ = L™, / e(^,x"l V^l ={x"',e}- (3.29) 

The left side of ()3.28p generates area-preserving diffeomorphisms with infinitesimal 
parameter e, where de is a harmonic 1-form on S. Thus we may write 



(3.30) 



where are time dependent functions such that de = iioj^ for some basis of 
harmonic 1-forms on S normalized with respect to the homology basis {6;}. The 
constraint (j3.28p is used in the quantum theory as a matching level condition. 



■^The constant is related to the total momentum (in the x direction) = /j, p_ by the 
relation = P^/vol(S). Here the volume is measured with respect to y/w . 



MEMBRANES IN LIGHT CONE GAUGE 



15 



We may now determine N and iV" so that the gauge conditions are preserved 
under evolution in t. Ensuring that = —ct by requiring that 



+ CT. 



(3.31) 



we find that 



In order that the second hght cone gauge condition be preserved, we require that 



N^^. (3.32) 



P-~cV^, /(A$ + A°$a) 



(3.33) 



which leads to the condition 

daiV^N""] = 0. (3.34) 

Let ri = ■^\/w £ab N°'d<j^; the condition ()3.34|) is equivalent to dVl — 0, i.e. that 
U, is closed. Consequently, on each contractible domain in E, we have that is 
(locally) exact: Q ^ df for some locally defined function /. The one form il may 
be globally decomposed, uniquely, into its harmonic and exact parts 

= /l + cfAoxact- (3.35) 

In any local coordinates, ft = {ha + 9aAcxact) do'"^- Thus 

TV- ^ ^ (hb + afcAexact) . (3.36) 



The form h may be expressed as h — dAharmonic for some multi-valued function 
Aharmonic; Write A = Aharmonic + Acxact- Notc that while A is multl-valued, dA is a 
well-defined geometric object. It is also useful to keep in mind that dx"^ are exact 
one forms, as they are single- valued functions E . This need not be true if 

the target space M.^ is replaced by a manifold with non-trivial topology. 

The light cone action may be obtained from (|3.22p by noticing that p.l9p - (|3.20p 
hold and thus 

§ = / Pf.drX'^= [ (pn^drX^^ - Cp+ + dr[p-X-]) , (3.37) 
JM JM 

where we may drop the last term as it is a total derivative. The reduced Hamiltio- 
nian in light cone gauge is then given by 

J{=^ (lPlI01 + l\{x,x}\A +p„{A,a;™}, (3.38) 



2 

where — {x™, x„} and we have made use of the the expression 

for 7. 



Notice that the constraints (|3.19p - (l3.20p are implemented in the action by the 
reduction procedure. The last term in is well-defined as it is expressed in terms 
of dA. In contrast, the expression A{p„i, a;™} is ill-defined. We now verify that the 
term pm{A,a;'"} indeed corresponds to a Lagrange multiplier term multiplied by 
the constraints. 
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The harmonic one form h may be expressed in terms of a basis {oj*} of harmonic 
l-formfl The basis contains 2g elements, where g is the genus of E. Thus 

h^\u\ i^l,...,2g, (3.39) 

where the coefficients \i are functions only of r. Making use of the bilinear Riemann 
identities (see, for example [2^) we have 

/ p„,{A,x™} = - / Ae,act|^,x™lv^+A, / ^dx^ (3.40) 

and arrive at the constraints (|3.27p - (|3.28p . We furthermore can interpret Acxact 
and the Xi as the Lagrange multipliers associated to these constraints; thus our 
remaining gauge freedom lies in the choice of these functions. 

We now turn to the equations of motion associated to the reduced Hamiltonian 
([OS]) . Defining Dr = dr + {■,A}, they are 

D^a;™ = drx"' + {x"\ A} = ^ (3.41) 



Dr[^)^^ + l^,A}^{{x^,x^},x-}. (3.42) 



Note that the constraints ()3.27|) - ()3.28p are preserved under this evolution. 

This system transforms covariantly under area-preserving diffeomorphisms pro- 
vided A transforms appropriately. We see that under area-preserving diffeomor- 
phisms generated by (|3.27p - (|3.28p . we have 

d5K^ddrC,+d{Q,K} = d{DrO, (3.43) 

where C is e + e as above. The harmonic part of dK transforms as 

6K = drh. (3.44) 

The exact part of (JA transforms as 

exact \ -^^-harmonic 

}, (3.45) 
where we have used that d{t + e, Aharmonic} is an exact 1-form. 



4. The initial value problem 



The gauge freedom present in the system allows us to fix Acxact and the Ai. 
Making the simple choice of setting each of these functions to zero, we see that 
Hamilton's equations of motion reduce to the second-order system (see also [i?! ]) 

dlx'"^ ^{{x^,x^},x^}. (4.1) 

It is interesting to note that under these choices, which fix and e up to time- 
independent parameters, the coordinate function t is in fact harmonic: D^t — 0, 
i.e., this the co- moving gauge. 

The remainder of this paper is devoted to studying the initial value problem 
for classical membranes in light cone gauge, as formulated in (|4.ip : For functions 



'Defined eg. with respect to w^b 
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{xo,uo) defined on E we show the existence of a function a; : [0,T] x E ^ M.^ ^ 
satisfying (14. 1|) and 

a;(0) = xo, drx{0) = uq. (4.2) 

We require that the initial data (a;o: ^^o) lie in an appropriate function space, as well 
as satisfy an appropriate version of the constraints in order that 

{drX"',X^} = (4.3) 

be satisfied by the corresponding solution x. 



4.1. The degenerate hyperbolic system (|4.ip . The main difficulty presented 
by the system (j4.1l) is that it is not strictly hyperbolic. First, note that for fixed x, 
the operator 



y^{{y"\x^},X^}^^d, 



w 



(4.4) 

V w L V 

is elliptic with symbol ^7(2;) 7(2;)°^. 

However, when the operator (|4.4p is applied to x itself, one must consider 

£ :.TK^{{x™,x"},x„}. (4.5) 

The first variation of £ is given by 

5L[x]y'^ - {{y™, x"}, x^} ~ {{y", x,,}, x^} + 2{{x™, x"}, y„}, (4.6) 

where we have made use of the Jacobi identity (|2.6p . The first term is the elliptic 
operator (|4.4p . but the second term also contributes to the symbol and thus bL^x) 
need not by strictly elliptic even if the metric 7(x)a& is Riemannian. 



Note that by the constraint (|4.3p . the second term in (|4.6p vanishes when we 
take y = drX. Thus the equation (|4.ip . together with the constraint (|4.3p . imply a 
non-degenerate hyperbolic equation for u — drX^ which we use below to construct 
the solution x. We also make use of the constraint when estimating solutions to 
the main equation (|4.ip . once their existence has been established. 



Returning to the operator £, we compute in local coordinates 

1 (4-7) 

= — (77°''(5mn - dadhx^ + lowcr order terms 

where e^„ —^ac^d Q^x^dd^n- Writing 

^^n^^l^T-'^mn-C'J (4-8) 

- \ (71^1^ - (e"" dax^vrf) (4.9) 



we have 



for vector VJ" e x R^~^. The first term appearing in L°^„ is diagonal and 
positive-definite when ^ah is Riemannian, but the second term can cause the symbol 
to be degenerate. For example, consider a (local) situation with x^ = ^ x^ = cP' ^ 
V2 = 1 and all other components of x and V zero. 
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These degeneracies associated to £ prevent us from constructing solutions by 
direct application of standard energy methods, as the energy-type quantities asso- 
ciated to L cannot be shown to adequately control approximate solutions. Thus we 
approach the initial value problem by considering a modified system, motivated by 
the linearization presented above, which we now describe. 



4.2. The modified system. Differentiating ()4.ip with respect to r, which has 
the effect of linearizing the system, and making use of the Jacobi identity ()2.6p we 
obtain 

dldrx^ = {{a.^E™, s"}, Xn} + 2{{a;™, rz;"}, drX^} - {{9.x", x„}, a;™}. (4.10) 

When the constraint (|4.3p is satisfied, the third term on the right vanishes; the 
remaining terms, when viewed as an operator acting on d^-x^ are non-degenerate. 
We take advantage of this structure in the following manner. 

For functions x,u let 

'D = 'D{x,u)={u^,x^}, (4.11) 

r ^ 2"\X,U) ^ drU"' - {{X^^ ,x''},Xn}, (4.12) 

={{w",a;"},x„} + 2{{a:"\x"},u„}. (4.13) 

A computation making use of the Jacobi identity and integration by parts shows 
that 



d 
dt 



J {1^^ + \3\^)V^ ^2j^(v{u"\drX^}+8m{dy'-A[x]drX"') 

+ {drX"',Xm}){drU"',drX„,})V^. (4.14) 



Thus if X, u satisfy 

l^d^u = A[x]u, 

then the conditions D = 0, 3™ = are preserved if initially satisfied. In particular, 
the function a; is a solution to the main equation (|4.ip satisfying the constraint 
(EJl). 



We thus approach the initial value problem for the equation (|4.ip by considering 
the modified system (|4.15p with initial data 

x(0)=xo, u{0) = uo, dru{0)^ui. (4.16) 

In order that the solution {x, u) to the modified system (|4.15p give rise to a solution 
of the main equation (|4.ip , we requir^ 

2)o = {<\xJfK„„ = 0, and (4.17) 
- - {{^0, ^0 }, <^}Sni = 0. (4.18) 

In the local existence result stated below, we require initial data with xq G i E 
^ and ui e H'^^^. The condition (|4.18p imposes an extra regularity condition 



■^One should also impose the global constraint condition (f^ uq ■ cLxq = on the data in order to 
have a solution to the Hamiltonian system; however, this condition does not play a role in our 
method for constructing solutions to the reduced equations lj4.1|l . 
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on xq. Note, however, that due to the degeneracy discussed in §4. II the condition 
(|4l8l) does not imply xq € H''+^. 

Provided ^{x)ab is non-degenerate, the operator A[x] is (quasi-diagonal) elliptic 
which in local coordinates can be written in divergence form: 

Xx) 



Aixy = da 



w 



■lixrdbu 



-ab 



^ab ^cd 
-cd ' 



dc 



dcx"^ddx"dtun 



hbddaU"' + 2daX'^dbX"ddUn) . 



(4.19) 



Note that by the anti-symmetry of the second line does not contain second 
derivatives of u. In particular, the symbol of A[x] is 'y{x) j{x)°''' . For notational 
convenience we denote the third line of (|4.19|) by b[x]u = b{x)°'daU, and the sum 
of the second and third lines by B[x]u — B{x)'^daU. Since j{x)ab — daX^^dbXn, the 
operator can be schematically written in two forms, 

da [77°''96u] +B[x]u - D [{Dx)^Du] 

+ {D^x){Dx)Du + [DxfDu 

= da [A'^^dbu] + h[x\u - D [{DxfDu] -f {DxfDu 

(recall D represents spatial derivative); here 



A 



ab 



A 



ab 



5,nnl{x) l{xf'' + 2 



^ab ^cd 



^c-^m^d^n 



(4.20) 
(4.21) 

(4.22) 



4.3. Well-posedness results. We show the existence of a solution {x,u) to the 
modified system (|4.15l) : by the discussion above if the constraints (|4.17p - (|4.18p are 
initially satisfied, this leads to a solution x to the main equation (|4.ip satisfying the 
constraint (14. 3p . This procedure implies no restriction on the solution, in particular 
all solutions of the membrane initial value problem with the appropriate regularity 
can be obtained in this way. 

Our approach views x = xq + Jq u as a. functional of u; thus x and u are required 
to have the same degree of spatial regularity. 

Theorem 4.1. Let fc > 4 and let xo,uo,ui G H'' x H'' x H''^^ such that ^{xo)ab 
is a non- degenerate Riemannian metric, and such that (|4.17p - (|4.18p are satisfied. 
Then we have the following. 

(1) There exists T > 0, depending continuously on the norm of the initial 
data and unique {x,u) € C^([0, T]; iJ'=) x C|. satisfying (|4J5l) - (|4J6)) . In 
particular, x G C^{[Q,t]; H^) is a solution to (|4.H) - ()4.2|) which satisfies 
gS]). 

(2) Let T* be the maximal time of existence for {x, u) as given above. Then 
either T^ — oo or 

sup (IIt^^IIl- + \\Dx\\w^.-- + ||I?u||l=o) = oo. (4.23) 

[0,T.) 
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Remark 4.1. An adaption of standard arguments, see [sj, §2.3] and citations therein, 
shows that the solution map iJ'^ xiJ'^'xi/'^^^ C^([0, T], i/*^) given by (ccq, uo, ui) 
(cc, u) G ^^([0, T], H'^), is continuous. Thus the initial value problem for the system 
(|4T7)) - (|4?T8ll is strongly well-posed. 

5. Proof of Theorem 14.11 

In this section we prove the existence of a solution (a;, u) as in point [1] of Theorem 
l4.1l to the initial value problem for the modified system ()4.15p - ()4.16p . and establish 
the continuation criterion stated in point Oof that theorem. First, we prove energy 
estimates for the linear system 

SI'S? - A[x]'^ = F, i<^,dr<S>)l^^^{uo,ui), (5.1) 

associated to (|4.15p . Once such energy estimates have been established, a standard 
sequence of arguments (see [s^l for example) implies the theorem. 

Below we use expressions like ||(9$o||l2, ||$o||; to denote norms calculated in 
terms of the initial data at r = 0. For the higher order norms, the higher order 
T-derivatives are calculated formally. 

5.1. The linear system. We consider first the linear system (j5.ip for some fixed 
X e C^([0, T]; iJ*-'), > 4, be such that the metric j{x)ab is non-degenerate on 
[0,T]. In this section we generally suppress the z-dependence of A, 7, and other 
quantities defined below. We furthermore denote by A any quantity which can be 
bounded by a constant times ||Da;||^;j^i,oo -|- ||_D97-a;||ioo . Let A^ = supjQ j,^ A. 

Define the energy E ^ E[^] 

E{t) - i^(|a,<i>p+7 7"'9a<i>™ab$™) V^. (5.2) 
Define A^; — Ke[x\ to be the smallest constant such that 

A^^lia^iii. <i?[*]<A|||a*||i. (5.3) 

on [0,r] for all v. Note that A^ is bounded when HDxHloo is. 

Differentiating with respect to r and integrating by parts we estimate 

a,s< lia,$iU2||F|U2 + A2p$||2,. (5.4) 

Integrating (|5.4p and applying Gronwall's lemma yields the standard basic energy 
estimate. 

Lemma 5.1. //$ G C^{[0,T]■L'^)C^C'^{[0,T]■,H^) solves dST]) with F e L\[0,T]; L"^) 
then 




This estimate is used below to construct estimates for ||<I>||; using the identity 
d^{D^dr^)- A{D^dr^) = D^drF+ [D^,A] dr^ + i[dr,A]<i>). (5.6) 



MEMBRANES IN LIGHT CONE GAUGE 



21 



Applying the energy estimate ()5.5p to the identity (j5.6|) yields the following, 

+ £ {\\drF\\H. + II [D\A] a.$|U. + \\D' i[dr,A] <&) lU.)) . (5.7) 

In order to obtain an energy estimate which closes, we need an estimate for 
spatial derivatives of This is accomplished by use of the elliptic estimate (|2.19p 
which implies 

WD'^L^ < A {\\D'-H\\l^ + \\AD'-H\\l2) . (5.8) 
Making use of the linear equation (|5.ip we have 

\\D'^\\l2<a{\\D'-H\\l2 + \\F\\hi-2 

+ \\D'-^d'ML^ + \\ [D'-\A]nL^) (5.9) 

Combining this with (|5.7p yields 

11*11/ < (1 + AT)A^e/o^A'(||<l>o||/ + ||i^(0)|U,-. + II [D'-^,A] $|U. 

+ / (||a,F||H,-2 + ||$||z_i 
Jo 

+ 11 [D'-\,A]drnL^ 

+ \\D'-'i[dr^A]<f)h^)) 

We estimate the commutator terms as follows. 
Lemma 5.2. Suppose ^ G , then 

\\[D,A]nL- < M\Dx\\H4nH2 

\\D{[dr,A]^)\\L'^ < A{{Dx))3\mH'^ 

If furthermore e for 3 < I < k, then by (^^71) Z?* G L°° and 

\\[D'-'',A]^\\l2 < AiA\\D-i'\\H.-2 + \\Dx\\h^-i\\D^\\l^) 
\\D'-^ {[dr,A]-^) IU2 < A{A\\D-^\\hi-i + {{Dx))i-i\\D-^\\l^). 

Proof. Direct computation shows 

\\[D,AMl^ < \\(D'A){Dn\L^ + \\{DA){D^n\L^ + \\{Db){Dn\L^- (5.12) 
Application of (|2.2ip . (|2.18p and (|2.17p implies (|5.11ap . Similar considerations 
imply (|5.11bp . To see (|5.11cp write 

[D^^^,A]'i' = [D^-^D, A] D-^ - {DA)D^-^D^ + D'"^ (^^^^r) - bD'-'^D^. (5.13) 
By the product estimate (|2.18p we have 

\\A\\h'-i <A\\Dx\\h'-i, and \\b\\H' < A\\Dx\\h'-2 (5.14) 
the second of which implies 

\\bDnw-- < A (A||D*||h^-2 + \\Dx\\Hi-2\m\\L-) ■ (5.15) 



(5.10) 



(5.11a) 
(5.11b) 

(5.11c) 
(5. lid) 
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The estimate ()5.11cp follows from commutator estimate ()2.22p . Finally, the product 
estimate ([2TT81) and Sobolev inequality ((2171) imply l!5.11djl . □ 



(5.16) 



As an immediate application of (|5.1ip we have 

II [D,A]dML- < ^Dx\\HA\drnH- 

Furthermore 

II [D,A]nL-<K\\Dx\\H'^ (^||$o||h2 +^ ||$||3^ 

\\[D''^,A]nL-<K^[\\D'^^\\Hi-2+ WDdMw--^ (5.17) 

+ A||i^x||^,-i (^\\D^o\\h-+J^ WDdrnH-^ 

where we have used the Sobolev inequality (|2.17p in the last line. 

Applying the commutator estimates above to (|5.10p . followed by application of 
Gronwall's lemma, yields the following energy inequalities. 

Lemma 5.3. A solution $ to the linear system ()5.ip can be estimated for t £ [0,T] 
by 

mU < Ab(1 + AT)e^^-(i+^-)^-«<^^»>^'-((l + A||i?x||H3)||<i>o||3 

+ \\m\\H^+ \\drF\\H.) (5.18) 

and, for k > 4, by 

\m\k < AEil + K^)e-^^i^+''M{{Dx))),_,,^ ^ K\\Dx\\H.-.m4k 

+ ||F(0)||h^-2+ j \\drF\\H.-2). (5.19) 

□ 

5.2. Application to non-linear system. In order to apply the above energy 
estimates to the non- linear system (|4.15p we fix /c > 4 and (xq, wq, ui) G x H'^ x 
iJ'^-i. For T,R>Q let 

Vlr^i^^Cli, : |||«||U,T<i?,«(0)-iio,a.«(0)=?/i}. (5.20) 

For any v e Vj^j, let x{v) — xo + /J" v. There exists ifo such that the restriction 
T < Kq/R implies Ae[x] < 2A£;[a;o]. We restrict to such values of T; thus 7[a;(t')]a6 
is uniformly elliptic for all w S Vj^rp. 

Define <i> = to be the solution to 

52$-^[x(t.)]$ = 0, ($,9,$)|^^„ = (7/o,"i). (5.21) 

The linear system (|5.2ip is hyperbolic and it follows from our assumptions that the 
coefficients are sufficiently regular that standard existence results apply, see [3|. 
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Estimating 

A||Dx||^.-i < WDxf^,^^ < 4 (1 + \\Dx4]J,^^ + T^\\\Dx\\\H>^-i^r) (5.22) 

the energy estimate (I5.19p . together with the Sobolev inequahty (|2.17p which pro- 
vides control of A, imphes the following. 

Lemma 5.4. There exists Rq > depending on xq,uo,Ui, such that for each 
R> Rq there exists Tr > such that for all T < the map v i-^ ^[v] takes V^rp 
to itself. □ 



Fixing some such R, we now show that for a possibly smaller T > 0, the map 
V i-^- is a contraction with respect to the C|, norm. 



Lemma 5.5. Let i?, be as given by Lemma \5.4\ For a possibly smaller value of 
Tn, we have for each T < Tr that 

mv'] - ^V^]\h,r < \\\v' -V^\h,r (5.23) 

for all v^.v'^ e Vj^rp. 



Proof Let $[1;*] and a:* = x{v'') for i = 1, 2. Then 

92($1 - $2) _ _ $2) ^ _ ^[2.2])$2^ (5 24) 

Schematically, 

dr [{A[x^] - A[x'^])<i>] = D [{drA{x^) - drA{x^)) D^] 

+ D[{A{x^)~A{x^))Ddr<^] 

+ {drh{x^) - drbix"^)) D$ 

+ {b{x^) - b{x^)) Ddr-^ 

Using the mean value theorem, along with product estimate (|2.18p and the Sobolev 
inequality (|2.17p the norm of each of the first two lines in (|5.25p is controlled by 
C(R){{D{x^ — x'^)))2- A direct estimate yields the same bound for the latter terms. 
Applying the energy estimate (|5.18p . one can choose T, depending on i?, such that 
ti 1-^ is a contraction □ 

A standard argument (see [4^,[33|) using that the solution map v $[1;] is 
bounded in and a contraction in Cj. yields part 1 of Theorem 14.11 



5.3. Continuation criterion. We now establish the second part of Theorem 14. II 
Consider the solution (x, u) to (|4.15p - (|4.16p and suppose [0,T*) is the maximal 
interval of existence. Let 

= sup (llDxIlvi/i,- + + Ab[x]) . (5.26) 

[0,r) 

Note that A^ is finite if and only if ||Z3a::||L°= and |||7(a;)|^"'^||L°c are. 
Using (|5.1ip and (|5.16p applied to (|5.10p we see for any r < that 

\\u\\k < (1 + Kr)^Ee^«^' ((1 + ^^)\\u4k + J\l + A^l)\\u\\k) ■ (5.27) 
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Thus Gronwall's lemma implies that 

|||w||U,r < (1 + A2)2Ai,e^(i+A'^)'A- ll^olU. (5.28) 

If Kt, < oo and < oo, then one may extend x,u,drU to [0,T*] such that their 
restriction to r = satisfies the hypotheses of the local existence theorem. As 
this contradicts the maximality of T* , we obtain the second part of Theorem 14.11 

6. Improved energy estimate 

Having established the existence of a solution x G C^{[0,T]; H'') to the main 
equation (j4.ip . we are able to establish an improved energy estimate by making use 
of the constraint (|4.3p . In particular, we obtain the following. 

Theorem 6.1. Let k > 4. For a solution x G {[0,T], H'') (g^l) with initial 
data (xq,uo), the maximal time of existence T* depends continuously on \\xQ\\fjk + 
||-uo||ff'=-i- 

Remark 6.1. Comparing with theorem 14. li we see that regularity requirement on 
the initial data is less. In particular, in theorem 14.11 it is required that xq,uq 
are in H'^, while in theorem 16.11 the regularity condition is of the type usually 
encountered for hyperbolic systems. Thus, from this point of view, when taking 
the area preserving constraint (|4.3p into account, the degenerate hyperbolic system 
(|4.ip behaves very much like a hyperbolic system. 

From the extension criterion, we know that x may be continued so long as 
At is finite. The Sobolev estimate (|2.17p implies that can be estimated by 
||i9a;o||^fc_i — WDxoW^k-i + ||uo||^fc-i- We now establish an estimate for this quan- 
tity by deriving an energy estimate for the main equation ()4.1|) itself. 

Retaining the notation above, we have 

A ^ A[x] = C{\\Dx\\w^,^ + \\Du\\lo^) (6.1) 

where we let the constant C increase (independent of x) as needed and A^ = 

SUP[o,r) A. 

Following the discussion in iJ4.1l we write (|4.ip as 

d^Xm - L^^ndadbx" = F„ (6.2) 

where F„i — Fm{Dx). In order to estimate derivatives of x we make use of the 
identity 

dl{D^x) - L{D^x) = D\F) - {d\l\ x (6.3) 

where L = L'^^dadb and we have dropped the m, n indices for notational conve- 
nience. 

Define the energy associated to — L by 

£[y] = £M(r) = 1 1^ + L,t.P^)ay'")(42/")) ■ (6.4) 

We are of course interested in the cases y = D^x for / = 0, . . . - 1. 
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From ()4.9p the non-degeneracy of 7(0:) implies there is a constant Ag = Ag [x; T] 
such that on [0, T] 

£M<A|||9y||i. and \\dy\\l. < Al^ J^{\dry\' + ^\Dyf^) . (6.5) 

Note that since Ag essentially controls the ellipticity of 'y{x)ab, it is effectively 
equivalent to Ag above. 

A priori, the energy £[y] does not necessarily control ||9j/||^2. Rather 

||9y|li.<A|(£[y]+ZM) (6.6) 

where 

Z[y]^^Jjy^\x,r,}^V^. (6.7) 

The following lemma shows that the "error" Z can be controlled for y ~ D^x, if x 
satisfies the constraint (|4.3p . 

Lemma 6.1. Since the solution x e C^([0, T]; i?'^) satisfies the constraint (|4.3p . 
t/ien /or < / < fc — 1 we /lawe 



<e" (^Z[i^'a;](0)+^"A2||I?'aa;||i2 
<e^ (q|axo||^,+rA^|||i?'9a:|||i.,,: 



(6.8) 



Proof. Applying _D' to the constraint yields 

{D'u"\x^} = - ^ (6.9) 
i=i 

from which we see that 

drZ[D^x]<Z[D^x]+y2 {D^^'u'^^D^x^Y (6.10) 

When j = 1 or Z we have 

/ {D'-'u"\D'x^f ^ <\\Ddx\\l^\\D'dx\\l^. (6.11) 
When 2 < j <l ~ 1 the Holder estimate implies 

/ {D^-^u"\D^XraY ^ <\\D^-Wu\\l^4D^'^D'^x\\l2, (6.12) 

JT, 

with 1 = 1^1 = From the interpolation estimate (|2.20p we have 

\\D^-Wu\\l2^ < C\\D^-^Du\\l^\\Du\\lz} (6.13) 

\\D^-^D^x\\l2, <C\\D^-^D^x\\^^J\\D^x\\l^. (6.14) 

Thus 

drZ[D^x] < Z[D^x]+ A'^\\D'dx\\l2. (6.15) 

Integrating and applying Gronwall's lemma yields (|6.8p . □ 
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The estimate ()6.8p for Z, together with ()6.6|) . impHes that when 

Te^ < AJ^At^ (6.16) 

then 

\\\dD'x\\\l2,<CAE I WdxoWli + snp E,[D' x]] . (6.17) 

This facihtates the construction of the foUowing energy estimate. 

Lemma 6.2. When (|6.16p holds the first derivatives dx of solution x G C^{[0,T]; H'') 
can be controlled in L°°{[Q,T]- H^) for Q <l < k - I by 

\\\dx\\\Hi,r<CAEefo^'\\dxo\\H'- (6.18) 

Proof. A straightforward computation using integration by parts yields 

drE-[y] < Wdrvh^dly - Ly\\L2 + A^\\Dy\\l,. (6.19) 

In order to apply this to y — D'x we must estimate the right side of (|6.3p . which 
we write schematically as 

D'F - [D',L] x = D'F- [D'D, L^^] Dx + D' {{DLtn){Dx)) . (6.20) 

Applying the commutator estimate (|2.22p and product estimate (|2.18l) we have 

\\D'F-[d\l]x\\l2<K^\\Dx\\hi. (6.21) 

Thus integrating (|6.19p and using (|6.17p we have 

\\\dx\\\l,^^ < CAl{\\dxo\\H' + £ A^DxWl^) (6-22) 

Applying Gronwall's lemma we obtain (|6.18p . □ 

We now show that there exists T > 0, depending only on |j9a;o||_f/3, such that 
the solution x to (|4.ip is defined on [0, T]. The key is to establish estimates for Al- 
and Ag, and to ensure that (|6.16p is satisfied. 

In order to estimate Ag note that there exists depending on A£[a::o], such that 
Ag [x] < 2A£ [xq] whenever 

T < KA^j}. (6.23) 

For such T, the condition (|6.16p follows from 

Te^ < CA£[xo]"^Ay2. (6.24) 

Making use of the Sobolev inequality (|2.17p the energy estimate (16.181) implies that 
when (|6.23p is satisfied 

At < CA£[xo]e^'^^||axo||ff.-i. (6.25) 

Since (|6.23p - (|6.24p - (|6.25p all hold with strict inequality for T = 0, we may choose 
T > 0, depending on ||9a;o||_fffc-i , such that they continue to hold on [0, T]. 
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